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Abstract 

Since it is not possible to predict when a Gamma-Ray Burst (GRB) will 
occur or when Active Galactic Nucleus (AGN) flaring activity starts, follow- 
up/monitoring ground telescopes must be located as uniformly as possible 
all over the world in order to collect data simultaneously with Fermi and 
Swift detections. However, there is a distinct gap in follow-up coverage of 
telescopes in the eastern U.S. region based on the operations of Swift. Moti- 
vated by this fact, we have constructed a 14" fully automated optical robotic 
telescope, Goddard Robotic Telescope (GRT), at the Goddard Geophysical 
and Astronomical Observatory. The aims of our robotic telescope are 1) 
to follow-up Swift/ Fermi GRBs and 2) to perform the coordinated optical 
observations of Fermi Large Area Telescope (LAT) AGN. Our telescope sys- 
tem consists of off-the-shelf hardware. With the focal reducer, we are able to 
match the field of view of Swift narrow instruments (20' x 20') . We started sci- 
entific observations in mid-November 2008 and GRT has been fully remotely 
operated since August 2009. The 3a upper limit in a 30-second exposure in 
the R filter is ~15.4 mag; however, we can reach to ~18 mag in a 600-second 
exposures. Due to the weather condition at the telescope site, our observing 
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efficiency is 30-40% on average. 
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1. Introduction 

The 7-ray emission which comes from extremely energetic 7-ray sources 
such as Gamma-Ray Bursts (GRBs) and Active Galactic Nuclei (AGNs) 
(especially blazars) is believed to be produced by accelerated particles in a 
relativistic jet with bulk Lorentz factors of ~100 for GRBs and ~10 for AGNs. 
A GRB let is formed when a massive star (> 10M Q ) collapses into a black 
hole QHl. On the other hand, an AGN jet is a continuous outflow from 
an active supermassive black hole (10 6 ~ 9 M Q ) [e.g. @|. Although the energy 
scale and bulk motion of jets differ by an order of magnitude between GRBs 
and AGNs, their radiation processes are expected to be similar. Therefore, 
understanding radiation processes in the context of shock physics and particle 
acceleration using both GRBs and AGNs will provide a deeper understanding 
of the the fundamental physics in these extreme environments. 

Recent observations of prompt GRB optical emission by ground/space 
robotic telescopes are providing key data for understanding the radiation 
mechanism of GRBs. The observations of GRB 050820A by RAPid Telescope 
for Optical Response (RAPTOR; (lj) and GRB 061121 by Swift UV/Optical 
Telescope (UVOT; |4J) indicate that there are at least two components in the 
prompt optical emission. One component is optical emission correlated with 
the prompt 7-ray emission. The other component is a smoothly rising and 
decaying component during the prompt 7-ray phase. The first component 
could be interpreted as emission from an internal shock because of the similar 
variability between the optical and 7-ray bands. The second component 
could be due to an external shock interacting with the inter-stellar matter. 
However, the Robotic Optical Transient Search Experiment (ROTSE) has 
shown that there are a couple of cases where the early optical emission does 
not correlate with the 7-ray emission (e.g. [(jj ). The extremely bright prompt 
optical emission from GRB 080319B observed by TORTORA and 'PI of the 
sky' challenges the standard picture of the GRB emission model |5|. 

Blazars form a sub-group of radio-loud AGNs and show an extreme vari- 
ability at all wavelengths [9|. The most accepted scenario is that a rotat- 
ing supermasive black hole surrounded by an accretion disk with an intense 
plasma jet closely aligned to the line of sight is responsible for the blazar 
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emission. However, fundamental understanding of the radiation process in 
blazars requires extensive monitoring campaigns at all wavelengths. In par- 
ticular, properties of variability (including flares) and spectra from simulta- 
neous data in various wavelengths provides key information for a physical 
understanding of blazars. 

2. Motivation and Telescope Specification 

Since it is not possible to predict when a GRB will occur or when AGN 
flaring activity starts, follow-up/monitoring ground telescopes must be lo- 
cated as uniformly as possible all over the world in order to collect data 
simultaneously with Fermi and Swift detections. Based on the operations of 
Swift, however, we notice a distinct gap in follow-up coverage in the eastern 
U.S. region (Figure [T]). This fact motivated us to construct a fully automated 
optical telescope at the Goddard Space Flight Center. If there is no GRB 
to observe, this telescope will perform extensive optical monitoring of Fermi 
Large Area Telescope (LAT) AGNs to identify the variability on the time 
scale of a day Q • 

The Goddard Robotic Telescope (GRT) system consists of a 14" Celestron 
Optical Telescope Assembly (OTA), an Astro-Physics 1200GTO mount, an 
Apogee U47 CCD camera, a JMI electronic focuser (EV2CM) with a PC con- 
troller (PCFC), and a Finger Lake Instrumentation color filter wheel (CFW- 
1-8) with standard Johnson U, B, V, R, I and Clear filters. The quantum 
efficiency of the camera and the transmission curves of each filter are shown 
in Figure [2j With the focal reducer (Celestron f/6.3), a 20' x 20' field of view 
has been achieved. The observatory dome is an Astro Haven 7ft clam-shell 
dome. Opening and closing of the dome is controlled by the serial interface. 
Figures |3] and 0] are pictures of the telescope and the dome. 

The dome, the weather station (Davis Instruments Vantage Pro2), the 
webcams (Logitech Quick Cam Chat), and the rain sensors are controlled by 
another Linux PC. The dome is directly connected by a serial port to the PC. 
The Linux-based weather station control software (wvieWLf) and the webcam 
driver (spca£0) are used to control the weather station and the webcams. 
The signal from rain sensors (Asuzac Inc.) is recorded by the I/O board 



1 http://www. wviewweather.com/ 
2 http://mxhaard. free.fr/index. html 
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Figure 1: Distribution of the GRB follow-up telescopes (information from GCN Circular 
of the Swift GRBs from September 2006 to March 2007). 



(ArduincH), and used to close the dome automatically in case of rain. 

The telescope site is the Goddard Geophysical and Astronomical Obser- 
vatory (GGAO) which is about 1.5 miles northeast from the Goddard main 
campus. 

3. Automation of the Observatory 

The whole telescope system, including the telescope mount, the focuser, 
the CCD camera and the filter wheel, is controlled by the ACP Observatory 
Control Software^ An advantageous feature of the ACP software is to per- 



3 http://www.arduino.cc/ 
4 http://acp.dc3.com/index2.html 
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Figure 2: Left: The quantum efficiency of U47 CCD camera. The CCD chip is CCD47-10 
(back illuminated CCD) by e2V technologies. Right: The transmission curves of each 
filter. 
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form a short test exposure (five seconds in our case) after each slew and to do 
astrometry automatically using the Guide Star Catalogs! to find the center 
sky coordinate of an image. If the pointing is off, ACP requests the mount 
to re-slew so that the requested position will be exactly at the center of an 
image. Furthermore, ACP will automatically construct a pointing model by 
recording an offset between the requested position and the actual slew po- 
sition. This pointing model will be used for every slew to achieve a better 
pointing at the initial attempt. 

We developed a script to make an observation schedule for a single night. 
This scheduler calculates the elevation of all objects in the database. It 
selects objects which are visible from the site (an elevation has to be larger 
than 35° to avoid the light pollution), and then schedules objects at the time 
when they reach the highest elevation. The output of the script is a text file 
with lines of ACP commands. ACP commands consist of an object name, 
a set of coordinates for each object, filters, exposure times of each filter, a 
bin size of images and the starting time of the observation. Since it is a 
text file, it is easy to modify the schedule by hand if we need to schedule 
an object which is not in the database. The sky flat, bias and dark frames 
are scheduled at the end of observations. The script runs once a day in the 
afternoon. This schedule file is automatically submitted through ACP by the 
observation control script which we describe next. As of September 2010, we 
have 58 Fermi LAT AGNs, 70 variable stars, 33 near-by galaxies to search 
for supernova, and 37 guest observers' objects in the database. 

A script to control the whole observatory has been developed. The script 
calculates the Sun's elevation every 10 minutes and sends commands to ACP 
or to the dome at specific Sun elevations to perform actions such as open- 
ing/closing the dome, changing the CCD temperature, submitting the sched- 
ule or re-submitting the schedule in case of a termination of on-going obser- 
vation for various reasons. 

The GRB position information is received via a Gamma-ray bursts Co- 
ordinates Network (GCN) socket connection. Once the GRB position has 
been received, our custom GCN receiver code will calculate the elevation of 
the GRB at the site, and if it is observable (an elevation of a GRB has to 
be higher than 20° and also the Sun's elevation has to be less than —20°), 
it will send commands to ACP to quit on-going observations and execute a 



5 A user can specify either the Guide Star Catalog 1.1 or the USNO A2.0 catalog. 
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GRB observation. We collect 200 frames of 5 sec, 80 frames of 30 sec, 120 
frames of 60 sec and 72 frames of 100 sec exposure of the GRB field in the 
time order given. Based on a real-time observation of several GRBs, GRT 
took 47 sec on-average to terminate on-going observations and complete a 
slew to the GRB location. The initial exposure was made 66 sec on average 
after we received the GRB position information through the GCN. 

The automatic script for LAT AGNs has also been developed to do the 
photometry. The script performs the relative photometry using SextractoJ^l 
and PyRAF0 adopting three reference stars defined for each object. The 
aperture photometry (IRAF phot task) is obtained using a circular extrac- 
tion region of 5" for the source and an annulus with inner /outer radius of 
5.5"/6.25" for the background. The script also extracts the publically avail- 
able light curves of Fermi LAT@ Swift BAT@ and Swift XRT0 of the AGNs. 
Examples of multi-band AGN light curves produced by this script are shown 
in Figure |5j All of the AGN light curves observed by GRT are available to 
the public from the dedicated web page0 

In Figure we summarize the automation of the observatory as a flowchart. 



4. Performance 

We used the observations of GRB 100915 to estimate the detection limit 
of our telescope system. The first several hours had good weather conditions. 
The CCD temperature was — 20°C. The observations were performed in the 
airmass range from 1.13 to 1.27 (about four hours) collecting 423 frames 
in total with four different exposures (5 sec, 30 sec, 60 sec and 100 sec). 
All the light frames are in the R filter and corrected for the bias, the dark 
and the flat. First, we extract the sources in the image using Sextractor. 
Then, we perform photometry using the IRAF phot task. The aperture 
foreground radius is 1.6 times the averaged FWHM of extracted sources. 
The background annulus with inner and outer radius are, respectively, 1.75 
and 2.0 times the averaged FWHM of the extracted sources. A signal-to- 



6 http: / / terapix.iap.fr / rubrique.php?id_rubrique=91 / 
7 http://www. stsci.edu/resources/software_hardware/pyraf 
8 http://fermi. gsfc.nasa.gov/ssc/data/access/lat/mslJc/ 
9 http://swift. gsfc.nasa.gov/docs/swift/results/transients/ 
10 http://www. swift, psu.edu/monitoring/ 

11 http://asd. gsfc.nasa.gov/Takanori.Sakamoto/GRT/index. html 
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noise ratio of extracted sources is calculated using the IRAF phot result O 
The USNO Bl catalog Rl magnitude is used as the R magnitude of the 
extracted sources. The 3a detection limit has been defined by fitting the 
signal-to-noise versus R magnitude data, and finding the 3a crossing point 
to the fit (see Figure |7]). The 3a detection limit of a 30 sec exposure in R 
is ~15.4 mag. To understand the detection limit as a function of exposure 
time, we repeated the same analysis of the images with different exposures 
(5 sec, 30 sec, 60 sec, 120 sec, 180 sec, 300 sec and 600 sec). We stacked 
30 sec images using swarrjjf] to create images with >30 sec exposure. The 
relationship between the 3a detection limit and exposure time is shown in 
Figure [SJ We reached ~18 mag in a 600 sec exposure. 

In Figure IHJ we summarize the fraction of good observations for each 
month for four AGN observations. Although we experienced few months 
with >60% of good weather night, the typical fraction of observable nights 
is 30-40%. 



5. Current Status 

We started fully automatic operation in mid- November 2008. Although, 
at the beginning, we had conservative operations by having a duty operator 
every night, GRT is completely human free as of August 2009. We are 
currently observing ~85 objects and collecting ~350 images per night. 

Our optical contribution to the multi-wavelength observing campaign of 
the 7-ray outburst of 3C 454.3 has been published [2|. Further collaborations 
with the Fermi/LAT team for Mrk 421, Mrk 501, 1ES 1959+650 and MG4 
J200 11 2+4352 are on-going. 

We have performed 13 follow-up observations of GRBs (12 Swift GRBs 
and one INTEGRAL GRB). Three of them were observed within 10 minutes 
after the GRB trigger. Although there is no detection of prompt optical 



12 The signal and the noise are calculated by 

signal = SUM - AREA x MSKY, 



and 



SUM - ARE Ax MSKY , „ 

noise = \ — + 1.23 ARE A x STDEV 2 . 
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emission or an afterglow so far, we confirmed that GRT can achieve a 17.5 
to 19 magnitude limit in the R filter, and also respond automatically and 
promptly to a GRB through the GCN socket connection. 
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Figure 5: Multi-wavelength light curves of Mrk 421 (top) and 3C 279 (bottom). From top 
to bottom for each source, the light curve in the 0.1-300 GeV band by Fermi/LAT, the 
15-50 keV band by Swift/BAT, the 0.5-10 keV band by Swift/XRT, and the optical R 
band by GRT. 
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Figure 6: Flowchart of GRT automation. 
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Figure 7: Relationship between signal-to-noise ratio and the Rl magnitude of the USNO 
Bl catalog for the image of the field of GRB 100915 at a 30 s exposure in the R filter. 
The horizontal dash-dotted line shows the 3cr detection limit. 
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Figure 8: Relationship between 3ct detection limit and exposure time. The best fit relation 
(dashed line) is 3a upper limit = 14.1754 x log(exposure time)+13.9666. 
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Figure 9: Fraction of good observations for each month for four AGN observations: Mrk 
501 (top left), Mrk 421 (top right), 3C273 (bottom left) and 1ES 1959+650 (bottom right). 
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